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Abstract 

We have developed a concrete quantum simulation scheme and experimentally simulated a pair- 
ing model on an NMR quantum computer. The design of our experiment includes choosing an 
appropriate initial state in order to make our scheme scalable in near future, and the accomplish- 
ment of our experiment makes use of twice Fourier transforms so that our method is applicable to 
other physical models. Our results show that the experimental simulation can give the spectrum 
of the simulated Hamiltonian. Consequently, the potential power of a quantum computer on the 
simulation of complex physical systems is verified. 
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I. INTRODUCTION 



Quantum computer (QC) can offer attractive ability in accelerating the computation. 

The most appealing feature from the classical computer, as Feynman |l| noted, is that QC 

can be used to simulate the physical behavior of a real system. Lloyd jij later confirmed 

this idea in a two-state array. Recently, Somaroo et al. j3| presented a general scheme for 

the quantum simulation with a physical system to simulate another. At present, among 

the physical realizations of QC, nuclear magnetic resonance (NMR) has shown the greatest 

achievements (for a review see e.g., jj, Q|) though it uses the nuclear spin ensemble |fllz|- 

With NMR technique, some interesting physical issues have been simulated successfully, 

ri 

including a four-level truncated quantum harmonic oscillator (31, a natural decoherence of 

n □ 

a two-spin system [8|], a two-qubit correlation function of the Fano- Anderson model |9[, 
a three-spin effective Hamiltonian Q], and the migration of excitation in an eight-state 



quantum system 



111- 



Quantum simulation utilizes a physical system, commonly referred to a QC, which is 
easier to control and measure, to simulate the dynamics of another one, which is rather a 
complex system, also called the simulated system. Quantum simulation is expected to solve 
the problems in quantum many-body theory which is dealt with more difficult by the classical 
methods. The simulated system we will discuss is the BCS model rJ which is of wide 
interest in condensed matter and gives a phenomenal explanation to superconductivity. In 
BCS model, one of the key problem is to find the e nerg y gap, i.e., the energy required at least 



o excite a quasiparticle from the Fermi surface pjf]. In fact, some approximate methods 
I5I [r| have been proposed to calculate the energy gap for a given BCS Hamiltonian. 
However, if we realize the quantum simulation, one can get the accurate result within a 
polynomial time and resources 

In this paper, we show that QC can be used to simulate the property of a superconductor, 
such as the BCS model. Based on a concrete scheme of quantum simulation on NMR 
QC which includes a method for choosing an appropriate working initial state and the 
usage of twice Fourier transforms Q], we experimentally obtained the spectrum of the 
BCS Hamiltonian, which coincided with the theoretic expectations. In Ref. authors 
had mentioned that the energy spectrum can be obtained by a fast Fourier transform after 
the measurement. In this paper, we deepen this point and show how to implement the 
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measurement of the spectrum on an NMR QC. 



II. SCHEME OF QUANTUM SIMULATION ON AN NMR QC 

For a complete quantum simulation procedure, there are three mainly qualitative steps 
- Hamiltonian mapping, the experimental simulation, and again mapping back to the sim- 
ulated system. We will obey these three steps to describe our quantum simulation scheme 
of the BCS model in the following. 

Firstly, the mapping from the BCS Hamiltonian to the Pauli operators. 

The reduced Hamiltonian of the BCS model is a pairing Hamiltonian ^,17, 2p| . 



N N 

Hbcs = h(%2 Tp( n ™ + n ~ m "> + V S c ln c -m c -i c i\, (!) 

m=l m,l=l 

where e m denotes the free electron energy from the Fermi surface; V the coupling coefficient 
where it is simplified as a constant 21 1; n± m = c^± m c± m the electron number operators; 
c m{ c m) the fermionic creation (annihilation) operator; m = 1,2, ••■ ,N represent all of 
relevant quantum numbers, and the electron pairs are labelled by the the quantum number 
m and —to, according to the Cooper pair situation where the paired electrons have equal 
energies but opposite momenta and spins: m = (k , t) and — m = (— k , |). For a typical 
metal superconductor, h~e m ~ 10~ 2 eV, hV ~ 10~ 6 eV j20l | . 

Our interested physical property of the BCS superconductor is the spectrum of the Hamil- 
tonian (JTJ), because an important parameter in superconductor — energy gap can be followed 
from it j22|. The energy gap in superconductor is the energy difference between the ground 
state of the element excitation to its first excited state. A lot of analytical and numerical 
methods (see 



23j for a review) have been developed to obtain the energy gap through the 
diagonalization of the BCS Hamiltonian. However, these methods are complicated and use 
some unproved approximations. In this paper, we show that the quantum simulation on a 
QC can solve it in a practical and direct way with reduced resources. 

In Ref. 111, authors had mapped the Hamiltonian (JTJ into the qubit space based on 
the isomorphic algebras of spin-fermion connection. The mapped Hamiltonian is also known 
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as the spin-analogy, 



N N 

H p = *E Y< + 2 E + ( 2 ) 

m=l l>m=l 

where a^, <r y , a 2 denote the Pauli operators. 
Secondly, experimental simulation. 

The details of the experimental simulation will be presented in the next section. Here we 
emphasize two problems. 

One is how to choose a working initial state for the quantum simulation. In order to 
get the spectrum of H p , one should select an appropriate initial state which can lead to a 
nonzero absorptive spectrum after the whole procedure of quantum simulation. Note that 
it is nontrivial to find such a state. In a recent work [l8| . we proposed a general method of 
the selection of an initial state for iV-qubit quantum simulation of BCS model. This method 
guarantees the existence of the absorptive peaks in the H p spectrum and the correspondence 
between the peaks and the H p enengy levels is clearly known. 

The other problem is how to obtain the spectrum of a Hamiltonian by quantum simula- 
tion on an NMR QC. We propose that the simulation scheme has to include twice Fourier 
transforms (fts) — the first ft is from temporal-domain to frequency- domain of the NMR 
free induction decay (FID), and the second ft is of the NMR amplitudes of a peak. 

To validate this point, let us briefly review the procedure of the NMR measure (the first 
ft). NMR internal Hamiltonian in a lab frame is [4j 



N N 

2 

NMR measure process under such an internal Hamiltonian 



N N 

flis = ^(E w W+ E « J <^- ( 3 ) 

i=l j>i=l _ 

ian is 4J], 



N 

rlab-L 



S NM r(uj) a ft[Tr(e- iH l^ Pfin e iH » £ af% (4) 



i=l 



where we label the evolution time during the NMR measure as t; Snmr{u) the NMR 
frequency-domain spectrum; pf in the state before the measurement, i.e., the state after 
a unitary evolution, which is, in our case, the evolution with H p ; o~f the project operator, 
o~ + = <j x + iOy, Tr the trace operation; ft the Fourier transform applied for the NMR 
measured signal which is the function of NMR measure time t. 
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Eq. (J3J) implies that after the NMR measure, the spectrum is just that of NMR Hamil- 
tonian self, but not the spectrum of the simulated Hamiltonian. Thus whatever the initial 
state or the evolution is, the spectrum is unchanged after the NMR measurement. From only 
NMR measure, we can not get the knowledge of the simulated Hamiltonian. However, the 
amplitude in NMR spectrum includes the information of the H p evolution time. If we select 
a series of discrete H p evolution time, and measure the corresponding amplitudes of the each 
spectrum, then apply another ft on these amplitudes, the result will be the spectrum of the 
Hamiltonian simulated. A more formal mathematic derivation for this procedure refers to 

From the above discussion, due to the characteristic of NMR measurement, twice Fourier 
transforms are required to get the spectrum of the simulated Hamiltonian when performing 
the quantum simulation. The first ft is the NMR measure applied for the NMR measure 
time, which obtains the NMR energy spectrum, and the second ft is applied for the simulated 
Hamiltonian evolution time, which obtains the spectrum of H p . 

Thirdly, mapping back to the simulated system. 

After the simulation, we are able to get the spectrum of H p . However, we have to remap 
it back into the BCS physical system and formulate the energy gap, which is the final 
aim of quantum simulation. In other words, we would like to know how the H p spectrum 
corresponds to the energy gap in BCS model. The derivation of the relationship between the 
eigenvalues of H p and the energy gap can be found in Ref. by the diagonalization of the 
spin-analogy Hamiltonian submatrices |24| and numerical calculations |2(j] . The energy gap 

;o states |01) and 1 10) which span a subspace with 



is a function of eigenvalues corresponding 
one spin-up states, i.e., one Cooper pair [17|. Thus, if we obtain the eigenvalues of the H p 
through the quantum simulation, the energy gap in BCS superconductor can be deduced. 



III. EXPERIMENTAL IMPLEMENTATION 



We performed the experiment on a Bruker AVANCE 400 MHz spectrumeter, keeping 
the temperature at 300K. The spin system is 13 C-labeled Chloroform (Cambridge Isotopes) 
dissolved in 6d acetone. Spins 13 C and 1 H are the two qubits and labeled as qubits 1 and 
2. The coupling constant is J = 214.9Hz. 

Our simulation can be described with 4 stages. 
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1) Prepare for the working initial state. 
The simplest initial state for a two-qubit system is 



with a pseudopure state based on spatial averaging 



mUp) ini = |00) + |01>. We began 
271 . |28[ , and then apply a Hadamard 



gate on the second qubit, say H 2 , to obtain the working initial state, 



|00) ^ |00) + |01). (5) 
2) Realize the unitary transformation. 

In this paper, we consider a 2-qubit system and its evolution with H p . The Hamiltonian 



is 



H P = \\ero\ + e 2 a 2 z + V{o\<? m + o\a% (6) 
and the corresponding propagator is 

JJ p (r) = e~ iHpT / h = e~ i ^ £ia * +E2a * +v ( a * <T * +a « a y^ T ^ ^ 

where we label the evolution time with H p as r. Note that we want to investigate whether 
the NMR control technique is adapted to be a testbed for the quantum simulation. In order 
to check the level of experimental accuracy to perform this simulation without any error 
in the unitary decomposition, we choose a special case of the Hamiltonian (JZj) - set the 
parameters E\ = e 2 = e, so as to the four terms in Eq. ((7|) are commuted with each other. 
Then Eq. ((7j) can be decomposed exactly into 

Actually, even if E\ ^ e 2 , there is no any important difference for our simulation procedure. 

Experimentally we set both of the spins at their resonance frequency respectively and use 
^-rotations to realize the first two terms in Eq. (JHJ), and then, use J-coupling evolutions to 
realize the last two terms. With the experimental one-qubit and two-qubit gates, Eq. (jHJ) 
can be expressed into 
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FIG. 1: Quantum network to realize the evolution of 2-qubit BCS Hamiltonian. Horizontal lines 
represent qubits 1 and 2. Time goes from left to right. The boxes marked by X, -X, Y, -Y denote 
the 7r/2 rotations along the respective directions, and (9{) z (i = 1,2) the rotation along z direction 
with the angle 0j. The combination of two dots and a vertical line denotes the J-coupling evolution, 
i.e., the Ising gate, namely e - l 2 nJa l a z T 3 . 



T 



e ~ i i cr l . e - i 5 7rJ,J z <J z' r3 



j21 1 ,12 



e *5 ff x . e *f^ 



(9) 



where r 3 is the J-coupling evolution time; r 1; r 2 are the z-pulse widths; the r.f. powers u 1 
and u> 2 can be specified by the 90 degree pulse widths of two spins respectively j^J. The 
z pulse e~ l ^°~' z9t on spin i (i = 1, 2) with the rotation angle 0j is implemented using the 
combination of x and y pulses (|)* — (0)* — (f )*_ x . and the rotation angle is 0, = cj^. The 
corresponding network to realize U (r) is shown in Fig. Q 
The parameters in Eq. Q and Eq. © satisfy 



er, 7rJr : ? = Vr. 



j x , l — UJ 2 T 2 = ET, 7TjT 3 = V T. (10) 

According to the values of the parameters in BCS Hamiltonian (J2J), we set e m /2iT ~ 10 4 
Hz, V/2ir ~ 1 Hz in the experimental simulation. For the H p evolution time r in Eq. (JBJ), 
we select 64 different values — the initial value is 0, the maximal value is s, and the 
increment is j- s. The NMR experimental evolution time {t±, T2, t^} can be calculated from 
the evolution time r using the Eq. (jlOjl . That means we repeat the experiments for 64 times 



TABLE I: Discrete sets of values of the evolution time used in experiment. t\ is the z-pulse width 
of qubit 1 ( 13 C), r 2 is the z-pulse width of qubit 2 ^H), and T3 is the J-coupling evolution time. 
Note the periods of tl, t 2 , T3 are ~ 40 (is, ~ 40 /is, ~ 10 ms, respectively. When the values exceed 
their respective periods, their values can be shortened. See text for details. 



Ti r 2 r 3 



initial values 

increment 17.6 [is 18.0 [is 1.4805 ms 

final values 19.7 [is 20.1 [is 0.2489 ms 

with the 64 sets of evolution time {n, r 2 , T3}. In experiment, the receiver phases are set 
with the same phase as that of the measurement of the first set {r 1; r 2 , r 3 }. We obtain the 
relative amplitudes of the finial states by integrating the peaks. 

Since the decoherence time 3] (T 2 relaxation time) of the two qubits 1 H and 13 C is no 
more than 3.3 s and 0.4 s, respectively, the whole evolution time should be controlled within 
the decoherence time. However, the values of evolution time in the latter sets {t\, t 2 , r 3 } 
are so long to exceed the decoherence time. To solve this problem, we experimentally use a 
technique to shorten the evolution time based on the fact that both the ^-rotation and the 
J-coupling evolution are the periodic rotations with the period 2ir. Therefore, the z-pulse 
width and the J-coupling evolution time longer than one period can be shortened less than 
one period. 

More specifically, the period of z-pulse is 2n radian, and since we use the hard pulses in 
heteronuclear system, the value of one period is ~ 40 /is. The period of J-coupling evolution 
is 2/ J. In our system, this value is ~ 10 ms. Note that compared with the J-coupling 
evolution time, the z-pulse width is so small and can be ignored. So the experiment time 
mainly depends on the J-coupling evolution time. Further note that there are two segments 
of the J-coupling evolution in our netwok (Fig. Q), the whole time of our experiment is about 
20 ms, which is well within the decoherence time. In our experiment, the z-pulse width and 
the J-coupling evolution time are listed in table. |U 

After the evolution, the final state is 



) fin = e- tHpT W« = e~ l£T [00> + cos(Vr) |01) - i sin(W) 1 10) . (11) 
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This state implies that there are two peaks appeared in NMR spectrum - one corresponds 
to the transition 1 00) <-» |10), and the other corresponds to 1 00) <->• 1 01 ) depended on the 
nonvanishing coefficients. 

3) The 1st ft. Since NMR spectrum of the eigenvalues is unchanged, the positions of 
the peaks include no useful information. What we need is to record the amplitudes of the 
peaks. For the initial state we select, there are two peaks in the final spectrum, as shown 
in Eq. (fTTj) — one is the peak of 13 C, another is the peak of l H. We measure the peaks of 
1 H, and record the real and imagery parts of the amplitudes corresponding to the different 
evolution time through the integration of the peaks. By substituting Eq. (fTTj) into Eq. (J3J), 
the amplitudes of the peak of 13 C and 1 H theoretically satisfied the equations 



Ampi = i sin(Vr)[cos(er) + i sin(er)], (12) 
Amp2 = cos(Vr)[cos(er) +isin(er)], (13) 

respectively. The measured amplitudes fit the Eq. (|T3|) . The real and imaginary parts of 
the experimentally measured amplitudes are shown in Fig. El 

4) The 2nd ft. Apply a discrete ft on the amplitudes, and obtain the spectrum of H p , 
shown in Fig. El The distance of the two peaks is 2 Hz from the figure. The value is same 
as the theoretical result which can trivially diagonalize the Hamiltonian (J2J) and calculate 
the energy difference corresponding to the eigenstate {|01),|10)}. The difference of the 
eigenvalues is just the information we want to know from the Hamiltonian H p . Since H p is 
the mapped one from the BCS Hamiltonian, we can remap the eigenvalues of H p to get the 
energy gap in BCS Hamiltonian The errors mainly arise as a result of the imperfect 

pulses and delays as well as the variability of the measurement. 

IV. CONCLUSION 

In this paper, a concrete experimental scheme of quantum simulation on NMR QC is 
presented. Two important components in our scheme are the selection of an appropriate 
working initial state and the usage of twice Fourier transforms. Based on this scheme, we 
simulate the evolution of the BCS Hamiltonian with an experimental technique to reduce the 
system evolution within the decoherence time. The simulation results show the BCS energy 
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FIG. 2: The real (a) and the imaginary parts (b) of the NMR measured amplitudes. The solid lines 
are fitted to the theoretical expectation. The stars are the experimental data. Note the horizontal 
axes is the evolution time r, from to 63/27T s. 

spectrum agrees with the theoretical value. This provides a practical method to obtain the 
characteristic of the complex physical system, and shows the potential power of the QC in 
quantum simulation. 

It is believed that quantum simulation is possible to solve the quantum dynamics. The 
ability to perform an efficient quantum simulation requires many degrees of freedom (perhaps 
10~100 qubits in the next generation QC (l2|)- In this paper, we are putting our research 
on those issues associated with problems that are difficult for quantum many-body theory 
to solve on classical computer. Although a 2-qubit quantum simulation of BCS model is 
shown in experiment, the fact that even small numbers of qubits leads to new insight into 
the mechanism of QC how to perform quantum simulation, and permit accurate prediction 
of the dynamics of complex quantum system if the technique developed future. 

Clearly, a number of challenges for the efficient simulation of physical system still remain 
at present jjj: Can we find the physical results in the simulated system after the algebras 
remapping? How can we realize a complex network of multi-qubit system with the enough 
accuracy? Can we find the a physical simulation problem that a QC can solve but not 



10 



35 




frequency (Hz) 



FIG. 3: The spectrum of the Hamiltonian H p obtained by the second Fourier transform. The solid 
line is fitted to the theoretical expectation. The distance between two peaks shows the difference 
of the eigenvalues of H p is 2 Hz, which is same as the theoretical calculation. 

the classical computer today? Nevertheless, we are sure that the concept presented here 
establishes a very first step to simulate a real physical system on QC. 
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Note added. -After this work was completed, we became aware of a similar work appeared 
in quant-ph/0410106 by C. Negrevergne, R. Somma, G. Ortiz, E. Knill, and R. Laflamme, 
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under the title liquid state NMR simulation of quantum many-body problem. 
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